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Abstract 

The paper develops the method for construction of the families 
of particular solutions to the nonlinear Partial Differential Equations 
(PDE) without relation to the complete integrability. Method is based 
on the specific link between algebraic matrix equations and PDE. 
Example of (2+2)-dimensional generalization of Burgers equation is 
given. 
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1 Introduction 



Analysis of nonlinear Partial Differential Equations (PDE) is severe prob- 
lem in mathematical physics. Many different methods have been developed 
for analytical investigation of nonlinear PDE during last decades: Inverse 
Scattering Problem [[l], ||, §, |], || Sato theory [§, ||, [10], [□]], Hirota 
bilinear method [12|, [13, II, [IH], Penleve method |fl6| , |I7, |18|], <9-problem 
19, 20. 21], with some generalizations [22, 23. 24, 25, 26], nonlinear d- 



problem [E7|, 



the algebraic matrix equation |3T], |32 



recent modification of the dressing method based on 

A wide class of PDE (so-called 
completely integrable systems) has been studied better then others. Never- 
theless, there are many methods which work in nonintegrable case as well: 

HB3I. 



12, 13, 14, 15, 16, 17, 18, 30 



We represent the method for construction of families of particular solu- 
tions to multidimensional nonlinear PDE without relation to the complete 
integrability. This method is based on general properties of linear algebraic 
matrix equations. Essentially we develop some ideas represented in the ref. f| 
and recently in the ref. [|TT|]. 

General algorithm is discussed in the Sec. 2. Sec. 3 deals with exam- 
ple of (2+2)-dimensional system of PDE together with family of particular 
solutions. By using an appropriate reduction, this system leads to Burgers 
equation. This example is simple demonstration of the method. Conclusions 
are given in the Sec. 4. 



2 General results 

The algorithm represented in this section is based on the fundamental prop- 
erties of linear matrix algebraic equation, which is written in the following 
form: 

W = $, (1) 

where \1/ = {ipij} is N x N matrix U = [u±, . . . , un] t , $ = [<f>i, • • • , 4>n] t ■ Let 
us recall these properties. 

1. If \l/ is nondegenerate matrix, i.e. 

det ^ ^ 0, (2) 
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then equation ([!]) has unique solution which can be written in the 
following form: 



U = ^- 1 $. (3) 
Only nondegenerate matrices ^ will be considered hereafter. 

2. (consequence of the previous property) If $ = and condition (|2|) is 
held, then the equation ([!]) has only the trivial solution 

U = 0. (4) 

3. (superposition principle) Consider the set of K matrix equations with 
the same matrix 

Wi = $i, i = l,...,K. (5) 

Then for any set of scalars bt (k = 1, . . . , K), function U = J2k=i ^kUk 
is solution of the following matrix equation 

K 

W = J2 b k$k. (6) 

fc=i 

4. (consequence of properties 2 and 3) If columns $j are linearly depen- 
dent, i.e there are scalars a/-, k = 1, . . . , K, such that 

K 

J2 a ^k = 0, (7) 

then 

A 



5>^ = 0. (8) 

i=i 

Note that analogous properties of linear integral equation have been used in 



the classical dressing method based on the <9-problem JT9|, ^0 . 

To relate the matrix equation (HD with the system of nonlinear PDE, 
one needs to introduce the set of additional variables x = (xi, . . . ,%) in 
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matrices $ and U, which are independent variables of succeeding system 
of nonlinear PDE. M is dimension of x-space. For this purpose we use the 
following system of linear differential equations 

M fc $ + iV*tf = k = l,...,P, (9) 

where and Nk are linear differential operators with matrix coefficients; 
Ak are N x N matrices, which are constant for the sake of simplicity The 
number of these equations, P, depends on situation. Equations (|) introduce 
variables x into the matrices \I> and $. Due to the eq.@, elements Ui of the 
column U are functions of variables x as well. 

Operators and Nk may not be arbitrary. They have to satisfy two 
conditions: 

1. Overdetermined system is compatible. 

2. Operators Mj~ and Nk have to provide existence of differential operator 
M. with non-constant matrix coefficients such that 

p 

M($-W) = J2( M ^ + N k^)U k + W = ^B(U) =0, (10) 

k=l 
P 

B(U) = J2 A ^t + V, 
fe=i 

where Uk and V are some matrix functions of elements and their 
derivatives (see eqs. ( ]2"5| - |2"ED as an example). 

Since det\l/ ^ 0, the last equation means: 

B(U)=0, (11) 

which represents the system of nonlinear PDE on functions itj. The structure 
of this system is defined by the equations @. 

One can see that the matrix equation (|l|) is generalization of the linear 
systems used in the refs.|| [TT| . In fact, let us introduce notations 

^i = ^f, (12) 
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and rewrite the matrix equation ([I]) in the form of system of N equations: 

N 

W = 5>^f, i = l,...,N. (13) 

3=1 

Two evident reductions are possible (these reductions will not be used in the 
succeeding sections). 

1. If 

Vf = fl£^ W , (i) =C + V l) , i,j = l,...,N, (14) 

then system ( jT5[) is system of iV ordinary differential equations on the func- 
tion ip® which is basic for the algorithm represented in the ref. [§]. 

2. Introduce multi-index (3 = (Pi, . . . , Pm), = ^kLiPk instead of 
index j (M is dimension of x-space), and use notations 

M 

*f = n^ (o =^ w , m\=i,...,N, as) 

= || a || = jv+l. (16) 

Then the system ( |T3"D can be written in the form 

N 

ga^ii) = upd^t/)®, i = l,...,Q, (17) 

11011=0 

where Q is the number of terms in the right hand side of the eq . fll7|) , a takes 
any value such that ||a|| = N + 1. Eq. (|17|) is basic in the ref. [pTf . 

Hereafter, we will use eq. fll3D instead of matrix equation (P. Superscript 
% runs the values from 1 to N, unless otherwise specified. 

3 Example: (2 +2) -dimensional generalization 
of Burgers equation 

In this section we consider (2+2)-dimensional equations with following no- 
tations for independent variables: 

x = x l7 y=x 2 , t = x 3 , t = x a . (18) 
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For convenience of construction, we separate the first term of the sum in the 
eq. (HI: 

N 

^ = Uli ,f + Y^u k ^f . (19) 

k=2 

Introduce two differential operators L\ and L 2 : 

Lt = h + 5 n u lx + 5 12 u ly , h = d t + d 2 x + jid x d y , (20) 
L 2 = h + S 21 u lx + 5 2 2Ui y , h = d T + dy + j 2 d x d y . 

Let us apply operator L\ to the eq. (|i~9|) and investigate the result 

- - 6 11 u lx ^ - 6 12 u ly <f>® + + u 1 M?)+ (21) 

u lx {2^f x + 7l # ,) + Ti^i^P .+ 

AT 

^((L^f+v^ = 

k=l 

Our purpose is to represent the equation ( |21| ) in the form of linear homoge- 
neous equation 

TV 

£*W£°=0, ( 22 ) 
i=i 

where are some expressions of and their derivatives. Owing to the 
property 4 (eqs.(^) and (||)), this would mean that Fi k = 0, & = 1, . . . , N. 
The simplest way to result in the eq. (|22"D is introduction of the system of 
differential equations (||) on the functions ipj and ) which looks as follows: 

I * 
* W - r (^iVwtf\) = E«^' (2 3 ) 

II k=2 

N 

<t> {i) -^\ = E"^' (2 4 ) 

iz fc=2 



6 



N N 



k=2 k=2 
N N 

4\ = E^-D^ 4\ = 52 a (^\ J>2, (26) 

k=2 k=2 
N N 

k^f = 52MP> = j>l, (27) 

fc=2 fc=2 
N N 

= 52mI\ = 52 ( 28 ) 



fc=2 fc=2 



Here parameters ay, ay, fy, (3ijk have to provide compatibility of the overde- 
termined system (p3|)-(p8|). Equation (|23| ) can be simplified if one eliminates 
by using the eq. (p4|) : 



^? y -G# x = 52a 2k ^\ G= S ^--, a 2k = ^(a lk -a lk ). (29) 



The system (^)-@ will be used below. 

Now the equation ([H]) takes the form (|22"D with 

-Pll = Wit + ^Ixo; + 71^1x1/ + SnUlUxx + 5l2«lMly + (30) 

(2/?l + /3 2 Jl)u 2x + fi\l\U% y = 0. 

Expressions for Fx/t with k > 1 will not be used. 
Analogously, apply operator L 2 to the eq.([l9|): 

- l 2 <f>® - 5nu u <j>® - 5 22 u ly ^ + (L 2Ul )^ ] + « 1 (i 2 ^f ) )+ (31) 

n l2/ (2^ i} y + 72 # J + 72«i,^V 

N 
k=l 

Due to the relations (|24|)-(f29"|) one gets analogous expression 

N 

52 F ^ = °> ( 32 ) 



3=1 
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if only 

72 1 

$n = -? (<*n7i - 2 <5i 2 ), $22 = ^(<5i27i72 + 2<Ju7i - 4<Ji 2 ). (33) 
7i 7i 

Coefficients F^k are expressed in terms of the functions Uk and their deriva- 
tives with 

72 

F 2 i = Ui T + Wi TO + ^Uixy + ^2(^1171 - 26i2)UiUi x + (34) 

^(5i27i72 + 25n7i - 46 12 ) 
7i 



Eqs.(|30|) and ( p4[) form the system of two nonlinear PDE on the functions u\ 
and it 2 . 



One shell discuss the structure of the system (24}-2~^). Operators L\ and 
L 2 , given by the formulas ( p0[) , determine the left hand sides of the equations 
in this system. The structure of the system of nonlinear PDE is defined by 
the operators Li and L 2 and by the right hand sides of the equations in the 
system (E3-E^). Only terms proportional to ipf* effect on the structure of the 



nonlinear system of PDE. At least one of the equations ( f2~4| - f2~9| ) must have 



this term (see eqs.(p5|)). Otherwise, only simplest solutions to the nonlinear 
system of PDE are available. Other terms enrich the family of available 
particular solutions. 

Function 0^ may not be linear combination of functions ij)^ (see eqs. 



(|23|j241) ). Otherwise all coefficients Uj in eq . (|T9|) are zero. Equations fl27|) 



and (p8|) introduce dependence on variables t and r. Right hand sides of 
these equations are zero in the examples below for the sake of simplicity. 

3.1 Construction of particular solutions 

This section is devoted to construction of some families of particular solutions 



to the system (^,gj). We will use simplified form of the system (|24| - P9| ): 

(35) 

(36) 

Mf + a, 2 ^\ (37) 



0«> 
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$12 


2 

7i 


4K 
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N N 

4\ = E a PH)rf. *? v = £«W)*V£°, 3>2, (38) 

k=2 k=2 

Z fc Vf = 0, Z fc (i) = 0, fe = 1,2, j = l,...,N (39) 

Below are two examples with N = 2 and N = 4. Example with N = 3 
does not have a principal difference with the case N = 2. Solutions from 
both families considered below are parameterized by the arbitrary functions 
of single variable. 



3.1.1 One-dimensional kink 

Let N = 2, superscript i takes the values 1 and 2. The system (|19j) has the 
following form: 

0« = u^? +u 2 ^\ (40) 



with 



ui = u 2 = — , (41) 



A 



4 2) ^ 



At 



A; 



4 2) 0( 2 ) 



(42) 



We need only equations (p5|)-(p7|) and ( p9|) (with j = 1, 2) to find functions 
ip^ (j = 1,2) and <f)^\ They admit the following solutions 



p = 


j cf(k)e kx+qy - 






f = 


j cf(k)e kx+qy - 






w = 


j c (i \k)e kx+qy - 







c w (A;) 



cfik) 
Pi 



(71/c 2 - a 32 7ifc + «i2/5i^i2), 
(A; - a 32 ), 



(43) 
(44) 
(45) 

(46) 

(47) 
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9 = £*, (48) 

k 2 

w = -k 2 -kq ll = - — {(3 1 +(3 2ll ), (49) 

Pi 

v = -q 2 -kq l2 = -^(f3 2 + f3 ll2 ), (50) 

Ml 

^ 12 $2 

Sn = t, — (#27i + 2/?i)> a 42 = -5-a 32 (51) 
Pili Pi 

One can see from the eqs. (|43| - |47|) that functions xp 2 are arbitrary functions 
of single variable (say, variable x). So that appropriate solutions ( f4"T|) depend 
on two arbitrary functions of single variable. 

As a simple example, let us take the following expressions for the functions 

ipf 1 and 

V>i = 3 (52) 

Pi 

m = fe^H S2£2 + hzpi^ (53) 

Pi Pi 

4 1] = s 1 E ll (54) 
# = s 2 E 2 + s 3 E 3 , (55) 



s E 

-Tl— (d!l2/5l5l2 - «327l^l +7l^l)> ( 56 ) 
P1A2 

S E 2 

(2) = -T^(ai2pl<5l2 - «327l^2 + 7l*l) + ( 57 ) 
Pl<5l2 



where 



^-^-{oLnPiSn - a 3 27i^3 + 71^3)) 
P1O12 

E = g* : ™ I +'"^+ w " i + l/ "' r n = 1 2 3 



q n = ^k n , Un = + M, ^ = -^(/? 2 + /?i72), (58) 

Pi Pi Pi 

n — 1,2, 3, 

which are in agreement with eqs.(K3|)-(BT]). Determinants are expressed 
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as follows: 

A = S -^((fa-k 2 )s 2 E 2 + (fa-k 3 )s 3 E 3 ), (59) 
Pi 

A 1 = -^^((A:! - A: 2 )(« 32 - ^ - fc 2 )s 2J E 2 + (60) 

Pl0 12 

(fa - k 3 )(a 32 -fa- k 3 )s 3 E 3 ), 

((fa - k 2 )(oL l2 i5 x b l2 - «3 2 7i + (x 32 ~i\(fa + k 2 ) - ^ifak 2 )s 2 E 2 + 
(fa - k 3 )(ai 2 (3i5i 2 - «3 2 7i + a 32 -ii(fa + k 3 ) - ^ifak 3 )s 3 E 3 ). 

Functions u\ and u 2 have no singularities provided that A 7^ for all values 
of parameters x, y, t, r, which happens if only 

sign((/ci - k 2 )s 2 ) = sign((fei - k 3 )s 3 ). (62) 



Taking into account relations (0), system ([3(].[l4|) gets the following form: 

uu + u lxx + -ixu lxy + - — (2/?! + p 2 'j 1 )u 1 u lx + (63) 

8l2UiU ly + (2/?i + f3 2 ^i)u 2x + Pxl-iVay = 0, 



720~12p2 _5l2 
71/ 

p2l2U 2x + (2fa + (3ll2)U 2y = 0, 



ui T + u lyy + i 2 u Xxy H —u x u\ x H 7r(2p2 + /3i7 2 )wiw la + (64) 

7iPi 7i Pi 



Note that owing to the relations ( ^8| - |50|) functions U\ and u 2 depend on 
two variables rather then four: 



X 



Pi 



T 



Pi 



-t + 



PI 



-T 



(65) 



Because of this fact, functions u\ and u 2} given by the formulas (|41|j59| - |6T 
represent single kink, which is essentially one-dimensional. Moreover, in 
virtue of eqs.(|65"D, the above system (B3LB4T) can be written in the form of dif- 



ferentiated Burghers equation on the function U\ with independent variables 
X and T. 
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3.1.2 Multidimensional kinks 



In this section N = 4, superscript i takes values from 1 to 4, unless otherwise 
specified. The system ( |i~9"D has the form 



(66) 



Formulas (^TJ) for the functions U\ and m 2 are held with determinants in the 
following form: 



A 



A, 



v4« 


if 


if 


if 


if 


if 


if 


if 


if 


if 


if 


if 


if 


if 


if 


if 


if 


0<» 


if 


if 


1 (2) 




if 


if 


if 




if 


if 


if 


4>f> 


ff 


if 



Ax 



4>m 


if 


if 


if 




if 


if 


if 




if 


if 


if 


4f.i 


if 


if 


if 



(67) 



Functions \ > ^ are defined by the same formulas (|35|-|3T|) and (R^ 



(j = 1,2) together with solutions ([ilf)- (|51~l) . Because of this, solution (f|l]) 
may depend on four arbitrary functions of single variable (compare with the 
paragraph below the eq.(|51~D). 

Functions r/>g and tp^' are solutions of the equations (|38| , |39D with j = 3, 4: 



',(0 



1>. 



(i) 

3 y 



Z fc Vf = 0, J = 3,4, fc 



1,2. 



"63^ + "64 + "62^2, ( 68 ) 

(69) 
(70) 



Q!83V'3 i) + "84^ + «82V'2 l) ) 



This is a system of linear nonhomogeneous equations with solution in the 
following general form: 



4 ] 



(71) 



where i/jjo an d V'S are general solution of homogeneous system, related with 



the system fl68H70| ), and particular solution of nonhomogeneous system (p8 - 
01) respectively. 
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The general solution of homogeneous system reads: 

,/,W _ „kix+qiy+u>it+i>iT _i_ „(*)»,( 2 ) k2X+q2y+&2t+i>2T ( 70\ 

Y30 — Pi '3 e +ft '3 e > 

■0® = pW r ( 1 ) e fcl^+9iy+^l*+l > l^ _|_ p(*) r ( 2 )gfe2a!+g2l/+W2t+52T ^73^ 

r W = (V^53>£ ) (74) 
"54 

<?j = («74&j + a 54 a 7 3 - a 53 a 7 4), j = 1,2, (75) 

"54 

/ a5 3 o;74 - 2a 5 3a 64 a74 + ag 4 a 74 + 4a| 4 a 8 3 



V"74 

where ujj and z/j are related with kj and by the equations 

Q i = ~ k ) ~ k^jli, Vj = -q) - kjW, (77) 



which follows from the eqs.([70|). 

The particular solution can be taken in the form (in virtue of the eq.([|4])) 

4p = Jr?(k)c^(k)exp((x + p 2 /p iy )k + uj(k)t + iy(k)T)dk, (78) 
4>® = J^ ) (k)c ( i\k)ex P ((x + p 2 /p 1 y)k + cu(k)t + u(k)r)dk, 

with u(k) and v{k) given by the eq. (^9||50D . Substitution of the eqs.(fTS|) into 
the eqs.( |68l|69|) results in the following expressions for Y\ and T 2 together 
with additional relations among parameters a^: 

ri(k) = — ^ , T 2 (k) = —Tx(k), (79) 

"52& - "52"53 - a 5A a 62 a 52 



"52 



a 5A a 62 (a 54: P 2 - a 74 p\ 



"52/32 a 62 (3 2 
"72 = —5 — , "82 = —5 — , (80) 
Pi Pi 

"74("84/3l - "74/^2) /olN 
"83 = 7 75 7JT2" X 81 

a 54 (o;54P2 - "74pl) 

(a 53 a U Pl - a 6i O!uPl + "54"84p\ - "53"54^2) 
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As an example, let us choose the following expressions for tpf^ and </>W 



4 4 



4 4 
(i 



4 4 



J. 1 ) „kiX+qiy+unt+Vi' 



1,2,3, 



S ^ C (A) e k4X+q 4 y+U4t+U4T _|_ s ^ c (5) e k 5 x+q 5 y+u} 5 t+v 5 T 
s gkiX+qty+Uit+ViT ^ 



1,2,3, 

k4X+q 4 y+uJ4t+U4T _|_ g^k^x+q^y+u^t+v^i 

= 1,2,3, 



«4e 

s , c (i) gkiX+qty+ujit+UiT ^ 



S4C 



(4) e k4X+q 4 y+uJ4t+V4T 
Jz\x+q\y+ib\t+v\T 



+ S 5 C 



(5) e k 5 x+q 5 y+u) 5 t+v 5 T 
^x+qxy+uijt+ujr 



ri(fc 3 )s 3 e fc3:c+933/+W3 * +iy3T , 

r^^e^ 94 ^ 4 *^ 47 " + rx(A;5)s 5 e fc5a;+95S/+OJ5 * +l/5T , 



(1) fcli+giy+iDit+i'i 



T + r 2 (/ci)sie 



r 2 (k 3 )s 3 e k3X+q3V+lJ3t+V3T , 

T 2 (k A )s 4 e k4X+q4y+UJ4t+V4T + T 2 (k 5 )s 5 e k5X+q5y+ul5t+U5T . 



(82) 
(83) 
(84) 
(85) 
(86) 
(87) 
(88) 
(89) 
(90) 
(91) 
(92) 
(93) 
(94) 
(95) 



where c^ n ' = c^ n \k n ), c x 

1 (n) 

1; 



, / 4 are given by the eq. (71) with r 3 



— Ci(k n ) are given by the eqs.(P|f|7) with = 
" 1 — 1; parameters k n , q n , ui n , v n are mu- 



tually related by the eqs. ([48| - |50"| ); relationship among parameters k n , q n , uj n , v n 
is given by the eqs. (|75l - [77|) ; T n are given by the eqs.(|79"D; s n and p n are arbi- 
trary parameters. 

General expressions for the functions u\ and u 2 are very complicated. As 
an example, let us fix the following list of parameters: 



h 

«12 

a 5 4 = "84 = 2, a 64 
Then, formulas (f|l]) read: 



1, k 2 = -1, k 3 ■ 
1, (5 2 = -3, 71 

"32 = «53 = «62 



2, /u4 — 2 ; 

72 = <$12 = 1, 
■ a 74 = 1, 



0. 



-1. 



ui = S( S5 e 2 ^ +12T )-2 S4 e 2(3 ^ 4i) ) x 



(96) 
(97) 
(98) 
(99) 



(100) 
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(343 S2 e^ + i* + 76p 2 e^ +i ^) , 
u 2 = t (l5680w^ + ^ + 8967 S2S5 e 2 ^ + f ^+i i + 24 -+ (101) 

4560 S4 P2e^ +%+8t+ ^ T + lMO^e^^) , 

A = 24010s 2 s 4 e^ + ^' + 21609s 2 s 5 e 2a:+ ^ + S' +24r + 
6840 S4 p 2 e^ +%+8t+ ^ T + 3420 S5 P2e^ x+1 ^ T . 

This solution has no singularities if 

sign(s 2 s 4 ) = sign(s 2 s 5 ) = sign(s 4 p 2 ) = sign(s 5 p 2 ) (102) 

The nonlinear system keeps the same form (^,^). The change of variables 
( |65| ) is not effective in this case, so that the nonlinear system is essentially 
(2+2)-dimensional. If one of parameters x, y, t or r is fixed, then functions 
Ui and -u 2 represent (2+l)-dimensional inelastic 3-kinks interaction. 



4 Conclusions 

The represented method is another way of using the algebraic system of equa- 
tions for analysis of nonlinear PDE. In comparison with method considered 
in 



30| , |31|) |32||, it supplies wide class of particular solutions for nonlinear PDE 
which may be solved by this method. In general, nonlinear PDE considered 
in this paper, are not completely integrable in classical sence. For instance, 
attempt to introduce commuting flows by the standard way leads, gener- 
ally speaking, to new constraints on the manifold of available solutions to 
nonlinear PDE. 

There is no formal restriction on the linear differential operators L\ and 
L 2 (see eqs.(^Uj)). This makes the algorithm fluxible, but its application field 
is not defined yet. An interesting problem is construction of commuting flows 
on submanifold of particular solutions to the given nonlinear PDE. 



References 

[1] C.S.Gardner, J.M.Greene, M.D.Kruskal and R.M.Miura, 
Phys.Rev.Lett., 19, 1095 (1967) 



15 



[2] V.E.Zakharov and A.B.Shabat, Funct.Anal.Appl. 8, 43 (1974) 

[3] V.E.Zakharov and A.B.Shabat, Funct.Anal.Appl. 13, 13 (1979) 

[4] V.E.Zakharov and S.V.Manakov, Theor.Math.Phys. 27, 283 (1976) 

[5] S.V.Manakov, Usp.Mat.Nauk 31, 245 (1976) (in Russian) 

[6] V.E.Zakharov, S.V.Manakov, S.P.Novikov and L.P.Pitaevsky, Theory of 
Solitons. The Inverse Problem Method, (Plenum Press, 1984) 

[7] M.J.Ablowitz and H.Segur, Solitons and Inverse Scattering Transform, 
(SIAM, Philadelphia, 1981) 

[8] Y.Ohta, J.Satsuma, D.Takahashi and T. Tokihiro, Progr. Theor.Phys. 
Suppl, No.94, p.210 (1988). 



[9] F.J.Plaza Martin, |math.AG/0008004 



[10] A.N.Parshin, Proc. Steklov Math.Inst., 224, 266 (1999) 
[11] A.I.Zenchuk, |arXiv:nlin.SI/0202053| vl (2002) 



[12] R.Hirota, in Lecture Notes in Mathematics 515, Springer- Verlag, New 
York (1976) 

[13] R.Hirota, J.Phys.Soc. Japan, 46, 312 (1979) 

[14] R.Hirota and J.Satsuma, Progr. Theoret.Phys. Suppl., 59, 64 (1976) 

[15] R.Hirota and J.Satsuma, J.Phys.Soc. Japan, 40, 891 (1976) 

[16] J.Weiss, M.Tabor and G.Carnevale, J.Math.Phys., 24, 522 (1983) 

[17] J.Weiss, J.Math.Phys., 24, 522 (1405) 

[18] PEstevez and P.Gordoa, J.Phys.A:Math.Gen., 23, 4831 (1990) 

[19] V.E.Zakharov and S.V.Manakov, Funct.Anal.Appl. 19, 11 (1985) 

[20] L.V.Bogdanov and S.V.Manakov, J.Phys.A:Math.Gen. 21, L537 (1988) 

16 



B.Konopelchenko, 5oteons in Multidimensions (World Scientific, Singa- 
pore, 1993) 

A.Degasperis, S.V.Manakov, A.I.Zenchuk, Physics Letters A, 249, 307 
(1998) 

A.I.Zenchuk, JETP Letters, 68, 715 ( 1998) 
A.I.Zenchuk, 119, 47 (1999) 

R.A. Kraenkel and A.I. Zenchuk, J.Phys.A:Math.Gen., 32, 4733 (1999) 

A. I.Zenchuk, J.Math.Phys., 41, 6248 (2000) 

B. G.Konopelchenko, L.Martinez Alonso and O.Ragnisco, 
|arXiv:nlin!si/0 1030231 vl (2001) 

L.V.Bogdanov, B.G.Konopelchenko, and L.Martinez Alonso, 

|arXiv:nlin.SI/0111062| vl (2001) 

B.G.Konopelchenko, L.Martinez Alonso and E.Medina, 
|arXiv:nlinlsi/0202013| vl (2002) 

A.I.Zenchuk, Phys.Lett.A, 277, (2000) 25 

A.I.Zenchuk, J.Math.Phys.,42, 5472 (2001) 

A.I.Zenchuk, J.Phys.A:Math.Gen. 35, 1791 (2002) 



17 



